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TIME-DOMAIN  ANALYSIS  OF  LUMPED/ DISTRIBUTED 
NETWORKS  FOR  EMC  APPLICATIONS 

1.0  INTRODUCTION 

Time-domain  analysis  of  electronic  networks,  which  deliberately  or 
inadvertently  incorporate  transmission  lines  and  other  distributed  elements, 
is  increasingly  important  in  EMC  applications.  Digital  systems  used  for 
control,  communications  and  computation  must  routinely  process  high  speed 
pulses.  Lightning,  EMP  and  various  other  effects  can  induce  unwanted  pulses 
into  both  analog  and  digital  systems.  Compatibility  analysis  requires  the 
ability  to  treat  pulses  in  such  systems  under  both  normal  and  abnormal 
conditions.  Shielding  requirements  and/or  potential  interference  levels  must 
be  evaluated.  Coupled  transmission  lines,  branched  cable  bundles,  and  other 
complex  combinations  of  wires  frequently  occur.  Transmission  lines  can  also 
be  used  to  model  a  variety  of  electromagnetic  shields. 

The  purpose  of  this  report  is  to  introduce  and  describe  a  new  technique 
[1]  suitable  for  the  time-domain  analysis  of  a  very  general  class  of  lumped, 
distributed  networks.  An  incidental  but  powerful  feature  of  the  new  technique 
is  that  the  same  system  equation  formulation  procedure  yields  either  time-domain 
or  frequency-domain  equations.  Considerable  saving  results  when  both 
frequency  and  time  solutions  are  to  be  obtained.  In  the  current  report  the 
basic  procedure  is  described  and  illustrated  with  examples.  The  analysis 
procedure  can  also  be  used  to  generate  time-domain  models  of  transmission 
lines  and  other  complex  structures.  This  feature  is  illustrated  by  generating 
an  exceptionally  simple  model  for  lossless  transmission  lines.  Finally,  a 

novel  concept  using  time-varying  reflection  coefficients  is  introduced. 

1 


It  should  be  emphasized  that  the  procedure  to  be  described  is  very 
general  and  applicable  to  a  broad  class  of  problems.  Further  applications  and 


in  particular  time-domain  models  for  lossy  transmission  lines  and  coupled 
line  structures  will  be  covered  in  future  reports. 


2.0  TIME-DOMAIN  ANALYSIS  OF  LUMPED/DISTRIBUTED  NETWORKS 

Time-domain  analysis  of  lumped  element  networks  is  well  established. 

Powerful  analytical  and  numerical  technqiues  are  readily  available,  including 
the  popular  state-space  and  Laplace-transf orm  methods.  General  purpose 
computer  programs  such  as  SCEPTRE  [2]  and  SPICE  [3]  provide  easy-to- implement 
time-domain  solutions  for  complex  lumped  systems  even  when  nonlinear,  time- 
varying,  and/or  active  elements  are  included. 

The  development  of  methods  for  transient  analysis  of  mixed  lumped-distributed- 
networks  is  of  relatively  recent  origin,  and  general  techniques  that  permit, 
for  example,  lossy  transmission  lines  of  arbitrary  lengths  and  nonlinear  active 
lumped  elements  are  not  yet  available.  Yet,  as  pointed  out  above,  the  time- 
domain  analysis  of  such  networks  is  increasingly  important  in  EMC  analysis  and 
prediction  as  well  as  normal  design  considerations  for  a  variety  of  pulse 
processing  systems  inlcuded  in  communications,  control,  computing  and  radar 
applications.  The  technique  to  be  described  is  widely  applicable  to  the 
analysis  of  such  systems. 

During  the  course  of  this  study,  a  substantial  literature  search  was 
carried  out.  The  most  pertinent  articles  and  books  are  listed  for  the  reader's 
convenience  [4-18].  While  the  technique  to  be  presented  is  significantly 
different  from  the  methods  found  in  the  literature,  the  present  concept  grew 
from  "wouldn't  it  be  nice  if  ••••"  considerations  following  a  reading  of 
Silverberg's  [4]  paper.  The  new  technique  has  been  successfully  applied  to  a 
wide  variety  of  problems.  The  impact  of  Silverberg's  work  is  gratefully 
acknowledged. 

2.1  SYSTEM  EQUATION  FORMULATION;  Part  I 

Assume  that  the  system  for  EMC  analysis  has  a  network  model  consisting  of 
interconnections  of  linear  distributed  elements,  dependent  sources,  and 

3 


T 


independent  sources.  Partition  the  network  into  two  parts  as  shown  in  Figure  1. 
One  part  consists  of  linear  (distributed  and/or  lumped)  elements.  The  other 
part  contains  any  lumped  nonlinear  or  time-varying  elements  and  independent 
sources. 


Figure  1.  Partitioned  network 

Silverberg's  [4]  procedure  is  to  solve  for  the  terminal  behavior  of  the 
linear  part  of  the  network  in  the  frequency  domain  and  then  convert  to  a 
terminal  time  domain  description  by  numerical  inverse  transform  techniques.  The 
time  domain  solution  for  the  whole  network  is  obtained  step  by  step  in  time 
at  the  interface  of  the  two  parts  by  simultaneous  solution  of  a  convolution 
equation  representing  the  linear  part  and  a  differential  equation  representing 
the  nonlinear  part.  The  simultaneous  solution  is  accomplished  at  each  time 
increment  by  solving  algebraic  equations  obtained  by  application  of  the  trape¬ 
zoidal  integration  rule  to  the  original  equations. 

For  the  moment  let  us  focus  our  attention  on  the  linear  part  of  the  network. 
Wouldn't  it  be  nice  if  the  frequency  domain  calculations  and  the  inverse 
transform  calculations  could  be  eliminated  and  all  calculations  be  performed 
directly  in  the  time  domain?  Computer  program  complexity,  memory  requirements, 
and  computational  time  could  all  be  significantly  reduced.  The  catch  is  that 
we  would  need  a  way  of  combining  element  descriptions  to  form  network 


descriptions  such  that  the  resulting  network  matrix  is  directly  compatible  with 
convolution  solutions.  Basically  this  implies  that  the  overall  system  matrix 
should  contain  only  sums  or  differences  of  individual  element  responses  (no 
products  or  quotients  allowed).  The  indefinite  admittance  matrix  [19]  appeared 
a  good  possibility,  but  because  of  the  type  of  systems  to  be  considered,  a  port- 
description  method  rather  than  a  terminal  description  method  was  desired. 

Kron's  transformation  methods  [20]  provided  the  inspiration  for  the  technique 
to  be  described.  Ultimately,  it  became  clear  that  for  the  present  requirements, 
the  formal  transformation  techniques  could  be  replaced  by  a  very  simple 
algorithm.  The  more  general  Kron’s  method  is  first  presented  and  then  the 
revised  and  simplified  algorithm  is  introduced  via  examples. 

At  this  point,  the  problem  statement  for  the  linear  part  of  the  network 
is  the  following.  Determine  a  scheme  for  representing  networks  such  that  given 
the  terminal  step-response  of  the  subnetworks  (or  elements)  in  the  time  domain, 
the  time-domain  terminal  step- response  for  the  connected  overall  network  can  be 
determined  as  sums  and  differences  of  the  individual  subnetwork  responses.  Then 
by  convolution  the  time-domain  terminal  response  of  the  overall  network  can  be 
determined  for  any  specified  set  of  inputs. 

2.1.1  Kron's  Method 

Short-circuit  admittance  parameters  will  be  used.  A  dual  impedance  formu¬ 
lation  has  also  been  used  successfully.  Following  Kron  [20],  let 

-*■  +-*■  ■+ 

I  =  YV 

be  the  given  matrix  equation  for  a  network,  where  I  is  the  current  vector, 

->■  •<— >■ 

V  is  the  voltage  vector,  and  Y  is  the  admittance  matrix.  Suppose  a  new  matrix 
equation  I  =  Y  V'  is  desired  for  the  given  network,  where  I  ,  V',  and  Y'  are 
the  new  current,  voltage  and  admittance  quantities,  respectfully.  Let  the 
relationship  between  the  old  and  new  voltage  quantities  be 


V  =  C  V' 


where  C  is  the  voltage  transformation  matrix.  In  most  cases  the  elements  of 
C  will  be  l's  and  0's. 

Power  in  the  network  must  be  the  same  for  either  choice  of  variables, 
since  the  network  is  in  no  way  changed  by  the  change  of  variables.  Thus, 


vtI  =  V*  I' 


must  be  true  where  subscript  "t"  indicates  "matrix  transpose"  and  superscript 
indicates  "complex  conjugate."  Substituting  V  =  CV'  yields 


V’ 

t 


C  I  = 
t 


V 

t 


'*p 


so  that 


and 


I* 


*-*■ 
=  C  I 
t 


+-*  <r-y  ■<-+ 

Y'  =  C  Y  C 


Collecting  the  above  results  yields  the  necessary  relationships  between  old 
and  new  network  quantities.  Associating  the  "old"  (unprimed)  quantities  with 
the  disconnected  subnetworks  and  the  "new"  (primed)  quantities  with  the  inter¬ 
connected  subnetworks  leads  to  an  algorithm  for  generating  the  system  equations. 
The  pertinent  equations  are  given  below. 

First  establish  the  relationship  between  the  "new"  and  "old”  voltage 

V 

quantities  to  generate  the  transformation  matrix,  C. 

V  =  C  V'  (I) 

Next  determine  the  "new"  admittance  matrix  from  the  "old"  admittance  matrix 
and  the  transformation  matrix  obtained  from  Equation  I. 

«-+•  «->*  e  ><  > 

Y'  =  Ct  Y  C  (II) 

The  "new"  system  equation  can  then  be  written. 
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The  "new"  currents  can  be  related  to  the  "old"  currents  and  the  transformation 


matrix. 

-*■  -*■*-*■ 

I'  -  ct  I  (IV) 

This  procedure  is  best  understood  through  examples. 

2.1.2  Combining  Subnetworks/Simplified  Method 

Initially  Kron's  method  is  used.  However,  it  will  be  shown  that  the 
constraint  of  allowing  only  addition  or  subtraction  of  subnetwork  matrix 
elements  leads  to  a  simple  algorithm  that  eliminate  a  number  of  steps  from 
Kron's  procedure.  The  underlying  feature  of  the  simplified  method  is  to  treat 
every  kind  of  connection  as  though  it  is  a  parallel  connection.  This  approach 
requires  that  open-circuited  ports  be  added  to  the  network  in  certain  situations. 
Such  additional  ports  are  like  ideal  voltmeter  connections  enabling  determination 
of  voltage  at  that  point  in  the  network  without  disturbing  the  system.  The 
added  open-circuited  ports  increase  the  size  of  the  system  matrix  but  the 
associated  currents  are  zero  and  the  overall  system  matrix  is  sparse.  The  net 
effect  of  this  type  of  transformation  on  computational  efficiency  has  thus  far 
seemed  to  be  increased  efficiency.  The  following  examples  illustrate  the 
conversion  from  Kron's  method  to  the  simplified  algorithm. 

2. 1.2.1  Example  1 

Given  the  2-port  networks  "A"  and  "B,"  each  represented  by  its  admittance 
matrix  as  shown  in  Figure  2a,  create  a  3-port  network  by  connecting  ports  1 
and  3  in  parallel  as  shown  in  Figure  2b.  Determine  the  admittance  matrix  of 
the  3-port  in  terms  of  the  original  unconnected  2-ports. 

Following  Kron's  method,  first  generate  the  "primitive"  admittance  matrix 
for  the  unconnected  subnetworks. 


Next  the  connection  matrix  relating  "old"  to  "new"  voltage  variables  is 
generated  as  in  Equation  I. 


(2) 


(3) 


Now  the  3-port  admittance  matrix  can  be  determined  using  Equation  II. 

-<-»-*  «->-  «— >■ 

Y 3-Port  =  Ct  YSUBC 
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(5) 


Notice  that  the  transformation  resulting  from  connecting  ports  1  and  3 
in  parallel  to  form  the  new  port  1*  yields  a  Y-matrix  which  could  have  been 

<  f 

obtained  by  adding  rows  1  and  3  and  columns  1  and  3  of  Y  to  form  the  1'  row 

SUB 
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and  column  of  Y.  „  .  Furthermore,  only  addition  of  individual  subnetwork 

3-Port 

■4 — ► 

elements  are  required  in  generating  Y3_port  f°r  connected  network. 

The  "new"  system  equation  can  now  be  written. 

V  2'  4’ 

1' 

=  2’ 

4' 

This  result  is  very  promising.  Is  there  some  way  other  types  of  connections 
could  be  looked  upon  as  parallel  connections  so  that  the  same  simple  results 
can  be  utilized  in  more  complex  situations?  To  test  this  idea,  consider  the 
cascade  connection. 


2. 1.2. 2  Example  2 

As  a  second  example,  consider  the  cascade  connection  of  the  two  2-ports 
as  shown  in  Figure  3a.  Common  practice  would  have  us  multiply  the  individual 
ABCD  parameters  to  obtain  the  new  ABCD  parameters  for  the  cascade  connection. 
However,  we  now  wish  to  use  Y-parameters,  to  avoid  products  and  quotients  of 
individual  terms  in  our  overall  description,  and  to  treat  the  connection  as 
though  it  is  "parallel"  if  possible.  This  can  be  accomplished  as  follows.  Add 
a  port  3'  in  parallel  with  port  3.  Notice  that  the  cascade  connection  shown 
in  Figure  3a  to  form  a  new  2-port  can  now  be  treated  as  a  parallel  combination 
of  ports  2  and  3  to  form  a  new  3-port  as  shown  in  Figure  3b.  If  port  3'  is 
open-circuited,  then  physically  the  networks  of  Figure  3a  and  3b  are  identical. 
However,  the  mathematical  descriptions  are  different.  In  the  first  case  the 
resulting  network  is  treated  as  a  2-port,  while  in  the  latter  case  it  is  treated 
as  a  3-port  with  1^  ■  0.  The  resulting  Y-raatrix  for  the  cascade  connection 
treated  as  a  constrained  (1^  =  0)  3-port  is  determined  as  follows.  First  form 
the  Y-raatrix  for  the  unconnected  subnetworks. 


0 

0 

0 

0 

vB 

vB 

Yn 

Y12 

vB 

yB 

Y21 

Y22 

Port  3  Open-circuited 

Figure  3. a)  Cascade  connection  of  2-ports 

b)  Cascade  connection  treated  as  a 
parallel  connection  with  added 
open-circuited  port 


The  rows  and  columns  of  Y  corresponding  to  ports  to  be  connected  in  parallel 

SUu 

are  now  added.  Ports  2  and  3  combine  to  form  port  3',  while  ports  1  and  4 
become  1'  and  4’,  respectively.  The  result  is  the  desired  Y-matrix  for  the 
cascade  combination  treated  as  a  constrained  3-port. 


The  more  conventional  2  by  2  matrix  representation  for  the  cascade 


u 


connection  can  be  obtained  by  eliminating  from  the  system  equations  (recalling 


R 


that  1^ 


0) .  The  resulting  Y-matrix  is 


(10) 


Which  obviously  includes  products  and  quotients  of  individual  2-port  terms, 
thereby  considerably  complicating  a  solution  by  convolution. 

2. 1.2. 3  Example  3 


True  parallel  connections  are  simple  and  require  no  added  open-circuited 
ports.  A  parallel  connection  of  one  port  of  a  3-port  network  with  one  port  of 
a  2-port  network  to  form  a  new  4-port  network  is  illustrated  in  Figure  4. 


2. 1.2. 4  Example  4 


A  series  interconnection  of  ports  in  terms  of  admittance  parameters  under 
the  constraint  that  only  sums  of  individual  subnetwork  admittance  parameters 
appear  in  the  result  requires  a  little  more  ingenuity.  An  auxiliary  connecting 
network  is  introduced.  The  series  connection  of  a  pair  of  ports  is  illustrated 
in  Figure  5  using  the  networks  of  Figure  4a.  Port  1  of  network  A  is  to  be 
series  connected  to  port  4  of  network  B.  A  "series  T"  is  connected  to  port  1 
and  the  Y-matrix  modified  as  shown.  This  operation  is  easily  done  automatically 
by  a  computer  upon  receiving  the  command  for  a  series  connection. 

2. 1.2. 5  Summary  Comments 

The  series,  parallel  and  cascade  connections  of  pairs  of  ports  permit  very 
general  networks  to  be  configured  from  subnetworks  (or  elements) .  The  very 
simple  procedure  outlined  in  these  examples  permits  system  equations  to  be 
formulated  treating  all  connections  as  though  they  are  parallel  with  the  result 
that  only  sums  of  subnetwork  admittance  terms  appear  as  desired. 

2.2  SYSTEM  EQUATION  FORMULATION:  PART  2 

Returning  now  to  the  total  network  consisting  of  linear  disbributed  and 
lumped  elements  plus  nonlinear  and  time-varying  lumped  elements.  The  network  is 
partitioned  as  shown  in  Figure  1.  The  solution  procedure  is  as  follows.  First, 
the  short-circuit  step-response  matrix  for  the  linear  part  of  the  system  is 
established  as  sums  of  the  individual  subnetwork  terms  as  described  in  the 
preceding  section,  then  a  matrix  convolution  equation  is  formed  relating  port 
voltages  and  currents  at  the  interface  between  the  linear  and  nonlinear  network 
parts.  The  interface  port  voltages  and  currents  are  simultaneously  constrained 
by  the  equations  for  the  nonlinear  part  of  the  network.  Both  convolution  and 
nonlinear  equations  are  represented  numerically  by  using  trapezoidal  (or  another 

appropriate  technique)  integration  leading  to  a  set  of  simultaneous  algebraic 
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equations  relating  port  voltages  and  currents  at  each  time  increment.  Solution 
of  these  equations  yields  desired  voltages  and  currents  at  each  time  increment. 

The  step-response  matrix  for  each  linear  network  may  be  determined  by 
measurement  or  calculation.  Let  A(s)  be  the  step-response  matrix  in  the 
Lap lace-transform  domain  and  Y(s)  be  the  short-circuit  admittance  matrix.  Then, 


A(s)  =  i  Y(s) 


r(t)  =  -1{A(s)  }  1{— Y(s)} 


(ID 

(12) 


where  a(t)  is  the  step-response  matrix  in  the  time  domain. 

Interface  port  currents  and  voltages  are  constrained  by  the  linear  part 

-V  -f 

of  the  network  as  follows.  I(s),  V(s)  are  vectors  of  port  currents  and 
voltages,  respectively. 


I(s)  =T(s)V(s) 

-  !Y(s>][s 

l(t)  =of_1{l(s)}  , 


s  V (s) 


(13) 


or 


i(t)  =  a(t)*v(t)U_1(t)  +  a(t)v(o  )  (14) 

where  v(t)  is  the  time  derivative  of  the  port  voltage  vector,  v(o  )  is  the 
initial  value  of  the  port  voltage  vector,  U  ^(t)  is  a  unit  step,  and  *  implies 
convolution.  The  nature  of  the  nonlinear  elements  is  assumed  to  be  such  that 
a  description  of  the  form 

v(t)  =  f(v(t),  i(t) ,  t)  (15) 

is  possible,  where  f(v(t),  i(t),  t)  is  a  matrix  whose  elements  are  explicit 

functions  of  v(t),  i(t)  and  t.  Equations  14  and  15  describe  the  network 

-► 

completely.  Given  the  initial  conditions  on  v(t),  we  can  in  principle  solve  for 
v(t)  and  i(t)  from  Equations  14  and  15.  Unless  the  matrix  of  functions  is 
extremely  simple,  the  solution  must  be  obtained  numerically.  Any  implicit 
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integration  technique  may  be  used.  Trapezoidal  (fixed  or  variable  interval) 
and  Gear  type  algorithms  [21]  have  proven  very  satisfactory.  For  ease  of 
presentation  the  fixed  interval  trapezoidal  method  will  be  presented. 

Let  A  be  the  interval  between  time  points.  Then  Equation  14  can  be 
written  as 

i(kA)  =  \  £ 
j  =  l 

xlv(jA)  -  v([j  -  1]A) I  +  a(kA)v(o+)  k  =  1,2, •••  (16) 


a( 


a([k  +  1  -  j]A)  +  a([k  -  j ]A) 


In  each  increment,  the  step  response  is  approximated  by  the  average  of  its  two 
end  point  values  and  the  derivative  of  the  voltage  is  approximated  by  the 
divided  difference  of  its  end  point  values.  Notice  that  i(kA)  in  Equation  16 

can  be  separated  into  two  parts,  one  depending  on  the  past  history  and  the 

->■ 

other  on  the  current  value  of  v(kA) ,  as  follows: 

i(kA)  =  i  (kA)  +  g  v(kA)  k  =  1,2, •••  (17) 

o  o 


where 


8o  2 


a(A) 


+1 

+  a(o^)j 


iQ(kA)  =  l  ^([k  +  1  -  j]A)  +*a(.[k  -  j]A)JJ^(jA)  -  v([j  -  1]A) 

-  VQv([k  -  1  ] A)  +  *a*(kA)v(o+)  (18) 

Thus, is  a  constant  matrix  equal  to  the  average  step  response  during  the 

first  time  interval.  The  vector  i  (kA)  can  be  treated  as  a  set  of  current 

o 

sources  whose  values  are  determined  by  the  past  history  of  v(kA).  For  a  given 

value  of  k,  i  (kA)  is  known.  In  effect,  a  lumped,  time-varying  terminal 
o 

equivalent  circuit  has  been  obtained  for  the  linear  (lumped-distributed)  part  of 
the  overall  network. 
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For  the  nonlinear  part,  from  Equation  18  we  have 


v(t) 


t 

* 

f(v(r),  i(i),  r)dx  +  v(t  -  A)  . 
t-A 


By  the  trapezoidal  integration  rule  we  have 
v(kA)  =  |  f (v(kA)  ,  i(kA)  ,  kA) 

+  |  f(v([k  -  1]A)  ,  i([k  -  1]A)  ,  (k  -  1)A) 


(19) 


+  v([k  -  1]A)  for  k  =  1,2, •••  (20) 

where  v(kA)  is  separated  into  two  parts,  one  depending  on  the  current  value 
and  the  other  on  the  past  history  of  v  and  i. 

The  solution  for  the  overall  network  is  obtained  by  solving  for  v(kA)  and 
i(kA)  simultaneously  from  Equations  16  and  20  at  each  time  increment  k  =  1,2,'*». 
Note  that  the  system  of  equations  is  algebraic  even  when  the  network  contains 
distributed  elements.  The  simplified  flow  chart  of  Figure  6  summarizes  the 
solution  procedure. 
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3.0  TIME-DOMAIN  EXAMPLES 

The  following  tutorial  examples  were  chosen  such  that  they  could  be 
verified  by  hand  calculations,  and  to  clearly  detail  the  solution  procedure. 


3.1  EXAMPLE  1 


Three  lossless  transmission  lines  are  interconnected  as  shown  in  Figure 
7a.  Determine  the  currents  i^,  i i^  given  that  v^(t)  =  U_^(t)  volts. 

The  network  is  first  redrawn,  breaking  the  circuit  into  subnetworks  whose 
step  responses  are  known  and  adding  open-circuited  ports  at  any  required  points 
as  shown  in  Figure  7b.  Let  A  =  1  usee.  This  information  is  supplied  to  the 
computer  causing  the  unconnected  subnetwork  matrix  to  be  established  as  given 
in  Equation  21  (zeroes  are  not  stored). 


aSUB 


(t) 


1 

4 

5 

2 

6 

3 

1 

A 

all 

A 

ai  -> 
12 

0 

0 

0 

0 

4 

A 

a21 

A 

a22 

0 

0 

0 

0 

5 

0 

0 

B 

311 

B 

a12 
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0 

6 

0 

0 

0 

0 

C 

all 

c 

a12 

3 

0 

0 

0 

0 

c 

a21 

C 

a22 

(21) 


The  a^'s  f°r  these  subnetworks  are  given  in  Figure  8.  A  variety  of  subnet¬ 
work  terras  frequently  needed  would  be  stored  and  available  in  a  general 


purpose  program.  New  a^'s  maY  input  as  equations,  tables  or  measured 
data. 


Next  the  interconnection  information  is  input  which,  in  this  case,  causes 
rows  and  columns  4,  5,  6  to  be  added  yielding  the  connected  network  matrix 
given  in  Equation  22.  The  individual  a^'s  are  given  in  Figure  8. 
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Figure  7, a)  Circuit  for  Example  1. 

Rq  =  50  fi,  T1  =  t2  =  T3  =  10  ^sec* 

v  (t)  =  U  (t)  volts 
8  1 

7.b)  Block  Diagram  of  Circuit  for  Example 
1  Showing  Added  Open-Circuited  Ports 
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(22) 


Notice  that  the  primes  have  been  dropped  from  the  port  designations  to  simplify 
writing  the  equations.  Interface  constraints  are  now  imposed.  In  general  this 
would  involve  a  set  of  equations  representing  the  nonlinear  and  source  part  of 
the  system.  In  this  example,  the  constraints  are  simply  v^  =  U  ^(t),  v^  =  0, 
v^  =  0,  i^  «  0.  Initial  conditions  are  v^o^)  =  1,  v2(o+)  =  v3(o+)  =  v^(o+)  =  0. 

Equation  16  can  now  be  written  for  this  example  as 


ix(kA) 

igiOcA) 

i2(kA) 

ig2(kA) 

i3(kA) 

103(kA) 

0 

Vk4> 

1 

2R 

0 

0 

0 

0 

1 

0 

1 

2R 

o 

0 

0 

0 

0 

0 

1_ 

R 

0 

0 

0 

0 

0 

0 

4 

R 

o 

v4(kA) 

(23) 


where  i^,  i^,  ig3»  ig^  must  be  calculated  at  each  time  increment  using 
Equation  18. 

For  k  =  1,  i^(A)  =  ig2(A)  =  i^^A)  =  ig4(A)  =  so  t^iat  from  Equation  23 
we  obtain 


24 


[  = 

il (A>  =  2R~ 
o 

i,(A)  =  0 

(24) 

i3(A)  =  0 
v4(A)  =  0 

No  change  occurs  in  any  variable  until  kA  =  T.  For  k  =  10,  i.e.,  kA  =  T, 
i01(10A)  =  iQ2(10A)  =  103(10A)  =  0,  and 
now  yields 

i^lOA) 

i2(10A) 
i3(10A) 
v4(10A) 

No  further  change  occurs  until  kA  =  2T, 

The  solution  proceeds  as  indicated  with 

3.2  EXAMPLE  2 

Let  the  network  of  Example  1  be  modified  to  include  a  nonlinear  element 

as  shown  in  Figure  10.  The  input  voltage  is  now  v  (t)  =  t  U  ^(t)  for  t  £  2  ysec 

and  v  (t)  =  0  for  t  >  2  ysec.  All  other  parameters  for  Example  1  remain 
S 

unchanged.  Determine  v2(t)  and  v3(t). 

The  setup  for  Example  1  remains  unchanged  except  for  the  new  input  voltage 
and  the  constraint  imposed  on  the  output  port  of  block  "C"  by  the  nonlinear 
device.  With  A  =  0.5  ysec  and  the  proper  nonlinear  constraint  imposed,  the 
program  yields  the  results  of  Figure  11,  which  can  be  easily  verified  by  hand 
calculations . 


i  . (10A)  =  -  ,  so  that  Equation  23 

U4  K 


_ i__ 

2R 

o 

=  0 

(25) 

=  0 

=  1 
=  4 

at  which  time  i^,  i2,  and  i3  all  change, 
the  final  results  shown  in  Figure  9. 
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Figure  10.  Circuit  for  Example  2. 

Rq  =  500  ft,  T1  =  t2  =  T3  =  10  ysec 

v  (t)  =  tU  , (t)  for  t  <  2  ysec 
g  -1  - 

=  0  for  t  >  2  ysec 


3.3  EXAMPLE  3 


A  network  consisting  of  three  sections  of  lossy  but  distortionless  trans¬ 
mission  lines  with  RC  loads  as  shown  in  Figure  12  is  driven  by  a  step  current 
generator.  Determine  v^(t)  and  V2(t).  This  example  is  one  used  by  Silverberg 
[4],  The  exact  solutions  for  this  problem  are 


v^(t)  =  (1  -  e  t)  volts 

V2(t)  =  (e  -  e  t)U_1(t  -  0.5)  volts  . 


(26) 


Results  computed  by  the  computer  program  are  shown  in  Figure  13  and  agree  with 
the  exact  results  to  six  decimal  places. 


Figure  12.  Circuit  for  Example  3. 

i  (t)  -  2U  ,(t)  Amps 
g  “I 

x  =  0.5  sec,  Yj_  -  Y2  “  Y3 


0  1.0  2.0  3,0 


Figure  13.  Solution  for  Example  3,  volt* 


4.0  FREQUENCY-DOMAIN  EXAMPLE 


As  pointed  out  in  Section  1.0,  the  same  procedure  can  be  used  to  formulate 
either  time-domain  or  frequency-domain  system  equations  which  utilize  only  sums 
or  differences  of  individual  subnetwork  matrix  elements.  The  resulting  system 
matrix  is  in  both  cases  usually  quite  sparse  with  the  associated  computational 
efficiency.  The  general  procedure  for  the  frequency  domain  is  illustrated  by 
the  following  example. 

Given  the  network  consisting  of  transmission  lines  connected  to  a  lumped 
element  junction  as  shown  in  Figure  14,  determine  the  average  power  dissipated 
in  the  40  ohm,  60  ohm  and  50  ohm  resistors.  Redraw  the  network  adding  open- 
circuited  ports  as  needed  (note  that  this  step  is  not  unique  and  a  number  of 
usable  combinations  can  be  generated).  The  result  is  shown  in  Figure  15. 

It  is  sometimes  convenient  (particularly  if  the  analysis  program  can  be 
used  on  an  interactive  basis)  to  create  the  final  system  matrix  as  a  sequence 
of  steps  involving  only  two  of  the  subnetworks  at  a  time.  This  procedure  is 
illustrated  with  the  current  example  as  follows. 

Given 
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to  be  combined  as  follows. 
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Zc  =  50  Q  j  '100  Q  j 


'  j  Lumped 
Junction 


Figure  14.  Circuit  for  frequency-domain  example. 

V  =  10  volts  peak, 
g 


so  n 

o-4V~°- 


Zc  =  50  n ijlOO  Q 


k-e*  =  ^-*1 


k-  B*  -  f  -4 


Figi.  e  15.  Block  diagram  with  added  open-circuited 
ports  for  circuit  of  Figure  14. 


From  Figure  15  the  output  ports  71  and  9'  are  short  circuited  and  ports  2',  4' 
and  5'  are  open  circuited,  so  that 


vy  =  =  0  and 


The  input  is  given  as  v  j  =  10  volts  peak.  The  system  equation  for  the  combined 
network  is  thus 


(33) 


These  equations  can  be  straightforwardly  solved  for  the  currents  ij,  i^,  and  i^ 

+->■ 

and  the  voltages  v'  v! ,  v'  since  all  elements  of  Y  are  known.  From  these 
2  4  5  new 

currents  calculate  the  power  delivered  to  the  resistors.  The  40  0  resistor 
receives  91.6  mwatts,  the  60  ft  resistor  87.9  mwatts,  and  the  50  ft  resistor  70.5 
mwatts . 
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5,0  TIME-DOMAIN  MODELS  FOR  TRANSMISSION  LINES 

The  techniques  of  Section  2.2  provide  not  only  a  very  general  analysis 
tool,  but  also  provide  a  powerful  scheme  for  generating  simple  time-domain 
models  for  complex  structures.  Modeling  will  be  illustrated  for  lossless 
transmission  lines.  Lossy  transmission  lines,  coupled  lines  and  other 
structures  will  be  covered  in  a  future  report.  These  models  may  be  used  to 
increase  the  power  of  existing  general  purpose  circuit  analysis  programs  such  as 
SCEPTRE  or  SPICE  by  providing  transmission  line  and  other  distributed  element 
capability  not  currenlty  available.  Alternatively,  the  models  may  form  the 
basis  of  new  analysis  programs. 

5.1  STEP-RESPONSE  OF  LOSSLESS  TRANSMISSION  LINES 

A  necessary  ingredient  for  time-domain  model  derivation  using  the  techniques 
of  Section  2.2  is  the  short-circuit  current  step  response  of  the  structure  to 
be  modeled.  This  step  response  may  be  determined  experimentally  or  analytically. 
For  the  lossless  transmission  line  the  step  response  may  be  determined  as 
follows. 

5.1.1  Unterminated  Line 

Figure  16  shows  schematically  a  lossless  transmission  line  of  length, 
characteristic  admittance  "G  and  propagation  factor  'jB."  Since  jB  =  jw£/v, 
where  U)  is  radian  frequency  and  v  is  phase  velocity,  jB?  corresponds  to  jBt, 
where  T  is  ?/v  equal  to  the  one-way  time  delay  through  the  line.  Next, 
transform  to  the  Laplace  transform  domain  replacing  jco  by  s  (the  Laplace 
variable) . 

The  admittance  matrix  for  the  lossless  transmission  line  in  the  Laplace 
domain  is 
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Y(s)  =  G 


coth  ST 

-1 

sinh  st 

-1 

sinh  st 

coth  s 

which  can  be  written  in  terms  of  exponentials  as  follows. 


Y(s)  =  G, 


ST  ,  -ST 
e  +  e 

-2 

ST  -ST 

e  -  e 

ST  -ST 

e  -  e 

-2 

ST  ,  -ST 
e  +  e 

ST  -ST 

e  -  e 

ST  -ST 

e  -  e 

(34) 


(35) 


Dividing  the  denominator  of  each  term  into  the  numerator  yields  the  following 
form. 


Y (s)  =  G 


— 

Yu 

Y12 

Y21 

Y22 

(36) 


where 

Y11  “  Y22  "  1  +  2  *e"2ST  +  '~4ST  + 

Y12  "  Y21  *  -2  {e"“T  +  e'3ST  +  «'5ST  +  "•> 

Taking  the  inverse  Laplace  transform  of  Y(s)  would  yield  the  time-domain  impulse 

1 

response.  The  step  response  is  obtained  from  the  inverse  transform  of  —  Y(s). 
Designate  the  step-response  matrix  by  a(t). 

^(t)  =  <3f-1  7  Y(s) 


an(t) 

a12(t) 

a2i(t> 

a22(t) 

(37) 


where 
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a..  =  a00  =  G  {U  _(t)  +  2U  . (t  -  2t)  +  2U  . (t  -  4t)  +  •••} 

11  11  o  -1  -1  -1 

a  „  =  a  =  -2G  {U  . (t  —  t)  +  U  . (t  -  3t)  +  U  (t  -  5t)  + 

11  11  o  -1  -1  -1 

and  U  (•)  Is  the  unit  step  function. 

5.1.2  "Properly"  Terminated  Lines 

If  a  "proper"  Gq  termination  is  used  on  the  input  or  output  of  the  lossless 
transmission  line,  the  step  response  reduces  to  a  simple  closed  form  rather 
than  the  infinite  series  of  the  unterminated  line.  The  properly  terminated 
form  has  some  advantages  when  used  with  SCEPTRE  and  SPICE.  If  a  termination 
different  from  G  is  desired,  an  appropriate  negative  or  positive  resistor 
is  inserted  in  series  with  the  Gq  termination.  Figure  17  gives  the  step 
response  for  input  and  output  terminated  lines. 

5.2  ASSOCIATED  MODEL 

Each  model  is  associated  with  (or  derived  using)  a  particular  numerical 
integration  algorithm.  Trapezoidal  integration  will  be  used  in  the  illus¬ 
trations. 

Consider  the  input  terminated  line  of  Figure  17a.  From  Section  2.2  (Eq.  14) 

i(t)  =  a(t)*v(t)  +*a(t)v(0+)  (38) 

For  fixed  time  increment  "A,"  the  discretized  form  of  Equation  38  is 

i(kA)  =  i  (kA)  +  g  v(kA)  (39) 

o  o 

~>  A — ► 

where  iQ(kA)  is  a  function  of  past  values  only  and  gQ  is  a  constant.  This 

— ►  — V 

means  that  only  i(kA)  and  v(kA)  are  functions  of  the  current  time  increment, 
kA.  The  discretized  equation  for  the  transmission  line  of  Figure  17a  takes 
the  following  form. 
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i02(kA) 
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(40) 


where 

Bu  ■  T  (,11(4)  +  *11<0+>1  '  T 

gl2-i[.l2(A)+a12(°+)].0 

*21  '  2  U21<A>  +  a21(0+>1  ’  ° 

822  "I  U22(M  +  a22(0+>1  ‘  Go 

Thus,  iQ1(kA)  and  i02(k4)maybe  interpreted  as  dependent  current  sources  whose 
values  may  be  calculated  from  results  obtained  at  prior  time  increments.  The 
simple  model  of  Figure  18  results. 


iL (A)  i2(kA) 


Figure  18.  Lossless  Transmission  Line  Model 

The  two  dependent  current  sources  can  be  evaluated  from  Equation  39.  For 
k  =  0,  v^(0)  and  v2(0)  must  be  supplied  as  initial  conditions.  For  k  >_  1 
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r— 

G  v  (0) 

i01(kA)  =  ~2 - [1+  U-l(kA  "  2T)1  "  v2(0)GQU_1(kA  -  T) 

G  G  k-1 

-  -f  v1((k  -  1)A)  +  -f  l  {[2  +  U_x((k  -  j  +  DA  -  2x) 

3=1 

+  U_J,((k  -  j)A  -  2t)  ]  [v^  (jA)  -  -  1)  A  )] 

-  2[U_1((k  -  j  +  DA  -  T)  U_LCCk  -  j)A  -  T)][v2(jA)  -  v2((j  -  DA)]}| 

(41) 

i  (kA)  =  -G  v  (0)U  (kA  -  T)  +  G  v  (0)  -  G  v,((k  -  1)A) 

1)2  O  1  -1  o  2  o  2 

G  k-1 

+  -f  l  (2[v  (jA)  -  v,((j  -  1)A)]  -  [U  .((k  -  j  +  1)A  -  t) 
j=l  Z 

+  ^((k  -  j)A  -  Dlfv^jA)  -  V;l((j  -  1)A)]}  (42) 

These  equations  are  easily  Implemented  yielding  an  efficient  transmission  line 
model  for  time-domain  calculations. 
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6.0  INSTANTANEOUS  REFELCT ION  COEFFICIENT 

A  very  general  concept  that  appears  quite  promising  in  dealing  with  complex 
circuits  is  that  of  a  time-varying  reflection  coefficient.  When  used  in 
conjunction  with  "fractional  time"  techniques,  the  time-varying  reflection 
coefficient  permits  certain  types  of  energy  relationships  to  be  deduced  in  a 
simple  fashion.  To  demonstrate  the  principle,  consider  the  circuit  shown  in 
Figure  19.  The  circuit  consists  of  a  lossless  transmission  line  driven  by 
a  triangular  voltage  pulse.  Generator  impedance  is  identical  to  the  character¬ 
istic  impedance  of  the  transmission  line.  Terminating  the  line  is  a  nonlinear 
resistor  whose  characteristics  (see  Figure  19b)  are  such  that  for  terminal 
voltage  amplitudes  less  than  V  the  resistor  looks  like  an  open  circuit.  For 
voltages  between  and  the  resistor  has  a  resistance  of  ohms,  i.e., 

equal  to  the  transmission  line  characteristic  impedance.  For  voltages  greater 
than  V2  the  terminating  resistor  looks  like  a  short  circuit.  Thus,  the 
reflection  coefficient,  p,  of  this  resistor  as  viewed  from  the  transmission 
line  will  be  a  function  of  the  terminal  voltage  and,  since  that  voltage  is 
time-varying,  the  reflection  coefficient  is  itself  time  varying.  In  fact,  the 
reflection  coefficient  is  p  =  +1  for  voltages  less  than  V^,  p  =  0  for  voltages 

between  and  V^,  and  p  =  -1  for  voltages  greater  than  V2«  The  fraction  of 

2 

incident  instantaneous  power  absorbed  by  the  nonlinear  resistor  is  1  -  p  . 

The  voltage  interval  -+  V2  during  which  power  is  absorbed  and  the  correspond¬ 
ing  time  intervals  for  a  triangular  input  pulse  are  shown  in  Figures  19c  and 
19d.  As  can  be  seen  from  the  sketches,  the  fractional  time  for  which  p  =  0 
decreases  as  the  input  pulse  amplitude  increases.  Since  in  this  case  power 
is  absorbed  only  when  p  =  0,  the  total  energy  (product  of  time  and  power) 
absorbed  by  the  nonlinear  resistor  decreases  as  the  input  pulse  amplitude 
increases.  The  energy  reflected  by  the  nonlinear  termination  is  absorbed  by 


43 


R 

ohms) 


Circuit 


0  V, 


-  - 


Stepwise 

Nonlinear 

Resistor 


Short  Circuit 


k— ■  At 


-H^tje- 


*iatk- 


Figure  19.  Instantaneous  Reflection  Coefficient  and 

Fractional  Time  Calculations  for  Predicting 
Bounds  on  Absorbed  Energy  (Explanation  in  text) 
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the  matched  load  at  the  source  end  of  the  line.  The  interesting  point  is 
that  for  this  and  related  type  v  -  i  characteristic  nonlinear  devices,  one 
can  say  that  if  the  input  energy  increases  by  a  factor  K,  the  energy  absorbed 
by  the  device  will  increase  by  a  factor  less  than  K,  i.e,  linear  extrapolation 
provides  an  upper  bound  estimate  for  energy  absorbed. 

Many  semiconductor  elements  have  v  -  i  characteristics  that  yield 
results  similar  to  that  of  the  above  example.  As  another  simple  example, 
consider  a  square-law  diode  on  the  end  of  a  transmission  line,  as  shown  in 
Figure  20a.  Since  the  v  -  i  relation  for  the  diode  is 

li  •  kvd  <“> 

where  k  is  a  constant,  the  instantaneous  admittance  of  the  diode  is 


y ,  =  —  =  kv ,  . 
J  d  v ,  d 


Diode  reflection  coefficient  referenced  to  the  transmission  line  is 

Y  -  kV 

*  T-TvT  <«> 

o  d 

and  after  some  algebra 

-  4kY  VJ 

1  _  p2  =  - (46) 

(kVd  +  V 

2 

A  plot  of  1  -  Pd  versus  is  given  in  Figure  20b.  The  percentage  of  incident 
power  absorbed  by  the  diode  first  increases  and  then  decreases  with  increasing 
terminal  voltage.  Maximum  absorbed  power  occurs  at  the  "match  point"  where 
the  apparent  resistance  of  the  diode  is  equal  to  the  characteristic  impedance 
of  the  transmission  line. 

Interpretation  of  power  absorbed  in  terms  of  time-varying  instantaneous 
reflection  coefficient  also  applies  to  circuits  with  reactance.  To  illustrate 
this  point  of  view  a  sequence  of  parallel  connected  ideal  lumped  elements  are 


considered  as  depicted  in  Figure  21.  Each  of  the  circuits  is  considered  to 
be  driven  by  a  sinusoidal  current  source  in  parallel  with  a  resistance.  This 
source  could  represent  a  Norton  equivalent  of  a  transmission  line. 

In  Figure  21a  a  single  resistor  is  driven  by  the  source.  The  instantaneous 
impedance  of  the  load  resistor  is  independent  of  time  and  its  reflection 
coefficient  is  a  real  constant.  Instantaneous  power  absorbed  is 


where 


2  2 

I  R  0 
,,  .  o  p  .  2 

PR  =  Vr  =  R  Sln 


RR 


(47) 


R  =  w 
p  R  +  R 


Time  average  power  absorbed  is 

7  2 

P  =  R  i  /2R 
avg  p  o 

Now  consider  what  happens  when  a  shunt  capacitor  is  introduced  as  shown  in 
Figure  21b. 

The  capacitor  alone  has  an  instantaneous  impedance  given  by 


(48) 


ZC(t)  = 


cv. 


V„  i  i  dt  V„ 


(49) 


where  V  =  time  derivative  of  V^,(t).  For  a  sinusoidal  input  the  instantaneous 

impedance  of  the  capacitor  varies  from  zero  (or  short  circuit)  to  infinite 

(or  open  circuit)  each  quarter  cycle  and  is  given  by 

1 

sin 


zc(t) 


u)C 


tan  u)t 


(50) 


Thus,  in  Figure  21b  the  constant  instantaneous  impedance  of  the  resistor  R  is 
now  paralleled  by  a  time-varying  instantaneous  resistance  that  varies  from  a 
short  circuit  to  an  open  circuit  each  quarter  cycle.  Since  the  time-varying 
reflection  coefficient  of  the  parallel  RC  combination  will  now  deviate  from 
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Figure  21.  Circuits  for  Demonstrating  "Instantaneous 

Impedance"  Conecpt  for  Sinusoidal  Excitation 
with  Reactive  Elements  Present 


that  of  the  resistor  alone  and  will  indeed  spend  a  significant  fraction  of  time 
dewelling  about  the  -1  value  corresponding  to  the  2ero  impedance  interval 


of  the  capacitor,  it  is  apparent  that  the  power  absorbed  by  the  resistor 
should  decrease.  Just  how  much  it  decreases  will  depend  upon  the  precise 
time-varying  instantaneous  impedance  of  the  capacitor  and  hence  upon  the 
value  of  capacitance  C  and  the  frequency  to  of  the  drive  signal.  Calculating 
the  instantaneous  power  absorbed  by  the  resistor  when  paralleled  by  a 
capcitor  yields 

2  2 
I  R 

pp(t)  =  - 2_E - sin2  (tot  -  toR  C)  (51) 

R[1  +  (toR  C)  ]  P 

P 

which  for  C  >  0  yields  a  decrease  in  amplitude  of  instantaneous  power  absorbed 
by  the  resistor.  Time  average  power  is 

P  =  R2I2/2R[1  +  (toR  C)2]  (52) 

avg  p  o  p 

which  again  shows  the  anticipated  decrease  in  power  absorbed  caused  by 
paralleling  the  resistor  with  the  instantaneous  impedance  of  the  capacitor. 

The  important  point  of  the  preceeding  discussion  is  our  ability  to  use 
time-varying  reflection  coefficient  to  predict  decreases  or  increases  in 
power  (or  energy)  absorbed  in  an  element  as  a  result  of  adding  some  new 
element  to  the  circuit  or  of  increasing  the  input  signal  level.  This 
technique  is  obviously  a  very  powerful  tool  where  bounds  or  limits  on  power 
(or  energy)  are  required. 

As  a  further  indication  of  the  power  of  this  "way  of  thinking"  about 
circuits,  consider  the  3  elements  in  parallel  as  shown  in  Figure  21c.  Here 
we  have  added  both  L  and  C  elements  in  shunt  with  the  resistor.  The  instan¬ 
taneous  impedances  of  both  the  capacitor  and  inductor  are  time-varying.  The 
formal  relationship  for  the  capacitor  was  given  above.  For  the  inductor  the 
instantaneous  impedance  is 


ZL(t)  =  L 


diL/dt 


(53) 


VLdt 


For  sinusoidal  input  signal  the  instantaneous  impedance  of  the  inductor 
varies  from  infinite  to  zero  each  quarter  cycle  and  is  given  by 


zL(t) 


=  iSL  cot  cat 
sin  ait 


(54) 


Thus,  in  Figure  21c  the  constant  instantaneous  impedance  of  the  resistor  R 
is  now  paralleled  by  a  pair  of  time-varying  instantaneous  impedances  one  of 
which  varies  from  open  to  short  and  the  other  from  short  to  open  over  each 
quarter  cycle  of  the  sinusoidal  signal.  The  net  result  is  that  the  instantaneous 
reflection  coefficient  of  the  combination  deviates  from  the  value  for  the 
resistor  alone,  and  for  every  case  except  one,  the  power  absorbed  by  the 
resistor  is  decreased  by  the  presence  of  L  and  C  in  parallel.  The  decrease 
in  power  absorbed  is  due  to  the  fact  that  the  instantaneous  reflection 
coefficient  of  the  combination  "dwells"  a  significant  fraction  of  time  about 
the  high  reflection  values  produced  by  the  reactive  elements  over  a  portion 
of  each  quarter  cycle.  Instantaneous  power  absorbed  by  the  resistor  is 


i2r2 

_°_E_ 


.  2 


R  ' 

2 

1  + 

wR  C  —  - J- 

_ 

wt  - 


r  R 

toR  C  -  — •*- 

P 


toL 


(55) 


which  is  less  than  the  power  absorbed  by  the  resistor  alone  for  all  cases 
except  the  case  where 


2 


w 


_1_ 

IX 


which  is  the  situation  we  normally  call  resonance.  From  the  point  of  view 
of  instantaneous  impedance  and  reflection  coefficient,  resonance  occurs  when 
two  time- varying  instantaneous  impedances  combine  in  such  fashion  as  to 
produce  a  resultant  constant  instantaneous  impedance.  Time  average  power 
absorbed  is 
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p 

avg 


R  I 
_JL°. 


f  R 

2 

1  + 

ojR  C - f 

p  0)L 

. 

(56) 


which  again  shows  the  anticipated  decrease  in  power  absorbed  (except  for  the 
resonant  case  wL  =  ^7  )  caused  by  the  shunting  effect  of  the  strongly  varying 
instantaneous  impedances  of  the  capacitor  and  inductor.  Note  clearly  that  it 
is  the  fraction  of  time  that  the  resultant  instantaneous  reflection  coefficient 
"dwells"  near  the  maximum  absorbing  point  that  determines  the  power  absorbed. 

The  above  discussion  and  examples  introduce  the  concepts  of  instantaneous 
impedance  and  reflection  coefficient  and  apply  them,  in  combination  with 
fractional  time  calculations ,  for  predictive  extrapolation  of  absorbed  power. 
The  concepts  are  exceptionally  general  and  may  be  applied  with  all  types  of 
elements,  linear  and  nonlinear,  active  and  passive.  While  it  requires  some 
time  to  become  accustomed  to  thinking  in  these  rather  unorthodox  terms,  the 
approach  seems  promising  for  general  purpose  analysis  of  a  number  of  EMC 
problems.  This  discussion  represents  only  the  beginning. 
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7.0  CONCLUSIONS 


A  new  technique  suitable  for  time-domain  analysis  of  a  very  general  class  of 
lumped/distributed  networks  is  introduced.  The  technique  is  useful  in  a  wide 
range  of  EMC  problems.  In  this  report  the  basic  procedure  is  described  and 
illustrated  with  examples.  Time-domain  models  of  transmission  lines  and 
other  structures  can  also  be  determined  using  the  analysis  technique.  Such 
models  are  useful  in  existing  CAD  programs  such  as  SCEPTRE  and  SPICE.  This 
feature  is  illustrated  by  generating  an  exceptionally  simple  model  for  lossless 
transmission  lines.  Finally,  a  novel  concept  using  time-varying  reflection 
coefficients  is  introduced. 
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